The Relativistic Dirac-Brueckner Approach to 
Nuclear Matter 
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An overview on the relativistic Dirac-Brueckner approach to the nuclear 
many-body problem is given. Different approximation schemes are discussed, 
with particular emphasis on the nuclear self-energy and the saturation mech- 
anism of nuclear matter. I will further discuss extensions of the standard 
approach, amongst other things the inclusion of non-nucleonic degrees of free- 
dom, many-body forces and finally compare relativistic and non-relativistic 
approaches. 

1 Introduction 

An ab initio description of dense nuclear matter which is based on QCD as the 
fundamental theory of strong interactions is presently not possible and will 
not be in the foreseeable future. The reason lies in the highly non-perturbative 
character of the formation of hadronic bound states and their interactions. 
Hence a quantitative description of nuclear many-particle systems has to be 
based on effective theories. Particularly successful are theories which's effec- 
tive degrees of freedom are hadrons, i.e. nucleons (and their excited states) 
and mesons. The nucleon-nucleon interaction is thereby described by the 
exchange of mesons as depicted in Fig. ^ Modern One-Boson- Exchange Po- 
tentials (OBEP), as e.g. the Bonn potentials PE], are usually based on the 
exchange of the six non-strange mesons: a (scalar, iso-scalar), uj (vector, iso- 
scalar), p (vector, iso-vector), tt (pseudo-scalar, iso-vector), rj (pseudoscalar, 
iso-scalar), S (scalar, iso-vector). 

The connection of the hadronic to the QCD world is reflected in the 
quark-hadron- duality. Due to the almost vanishing masses of the light (cur- 
rent) quarks chiral symmetry of the QCD Lagrangian is almost fulfilled. It is, 
however, spontaneously broken by the large non-vanishing vacuum expecta- 
tion values of the quark and gluon condensates which are responsible for the 
finite hadron masses. These condensates which are the basic quantities in the 
non-perturbative regime of QCD change dramatically in the medium. The 
connection to the energy density e in terms of hadronic degrees of freedom is 
in principle given via the Hellmann-Feynman theorem 
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where the sum runs over all hadronic contributions h to the energy density 
e. However, neither this part can uniquely be fixed nor the second part of 
eq. i.e. the derivatives of the hadron masses with respect to the current 
quark masses. Thus QCD can help to constrain hadronic theories but an 
exact mapping of the two worlds is still a dream for the future. 
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Fig. 1. Schematic representation of the One-Boson- Exchange model for the 
nucleon-nucleon interaction. The figure is taken from 

However, also within hadronic theories a perturbative approach to the 
strongly interacting nuclear systems is not possible. A systematic summation 
of diagrams up to infinite order in terms of the Brueckner hole-line expansion 
turned out to be an appropriate treatment. Already in lowest order which 
corresponds to standard Brueckner theory the saturation of nuclear matter 
can be described at least qualitatively [H IH] . In Brueckner theory the T- 
matrix (or Brueckner G-matrix) serves as an effective in-medium two-body 
interaction. It is determined by a self-consistent summation of the ladder 
diagrams in a quasi-potential approximation (Thompson equation) to the 
Bethe-Salpeter equation. The character of the bare nucleon-nucleon inter- 
action, in particular the repulsive short range part (hard-core) requires to 
account for two-body correlations in a self-consistent way. The effect of the 
correlations on the two-nucleon wave function in the medium is schematically 
depicted in Fig. [21 

However, non-relativistic Brueckner calculations are not able to meet the 
empirical saturation point of nuclear matter (psat = 0.16 fm~^, -Bbind = — 16 
MeV). In contrast, the saturation points obtained for various types of NN- 
potcntials were all located on the so-called Coester line 5 (for a recent review 
see j3j) in the e — g plane. A breakthrough was achieved when first relativistic 
(Dirac-) Brueckner-Hartree-Fock (DBHF) calculations were performed in the 
eighties |H1 El ■ Now the Coester line was shifted much closer towards the 
empirical area of saturation. One reason for the success of the relativistic 
approach is usually attributed to the fact that the dressing of the in-medium 
spinors introduces a density dependence to the interaction which is missing 
in the non-relativistic treatment. In the latter case the inclusion of three- 
body forces can lead to similar effects. The occurrence of many-body forces 
is, however, closely connected to the inclusion of non-nucleon degrees of free- 
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Fig. 2. Effect of the two-body correlations on the two-nucleon wave function as a 
function of the relative distance r. The nucleon-nucleon potential is schematically 
indicated. It shows the typical short-range repulsion { hard core) together with the 
intermediate and long-range attractive parts. The figure is taken from 



dom, i.e. resonances. A discussion of these aspects as well as a comparison of 
relativistic versus non-relativistic approaches will be given in Sec. 5. 

Relativistic Brueckner calculations are not straightforward and the ap- 
proaches of various groups [5] El d El El Ej are similar but differ 
in detail, depending on solution techniques and the particular approximations 
made. The intention of the present work is to review the standard relativistic 
Brueckner approach and to discuss implications of the several approximation 
schemes, as well as the role of higher order correlations, Pauli effects, and the 
special role of Dirac phenomenology. 

Large part of the present discussion will be devoted to the determination 
of the nuclear self-energy. To determine its Lorentz structure and momentum 
dependence, the T-matrix has to be decomposed into Lorentz components, 
i.e. scalar, vector, tensor, etc. contributions. This procedure is not free from 
ambiguities |16| . Due to identical matrix elements for positive energy states 
pseudo-scalar and pseudo-vector components cannot uniquely be disentan- 
gled for on-shell scattering. However, with a pseudo-scalar vertex the pion 
couples maximally to negative energy states which are not included in the 
standard Brueckner approach. This is inconsistent with the potentials used 
since OBEPs are usually based on the no-sea approximation. Hence, pseudo- 
scalar contributions due to the one-7r exchange (OPE) lead to large and 
spurious contributions from negative energy states. In jl2| it was shown that 
such spurious contributions dominate the momentum dependence of the nu- 
clear self-energy, and, in particular, lead to an artificially strong momentum 
dependence inside the Fermi sea. It was further demonstrated jT2] that pre- 
vious methods jSl ^2 used to cure this problem fail and finally a new and 
reliable method was proposed to remove those spurious contributions from 
the T-matrix jl3|. If calculations are performed in full Dirac space, i.e. in- 
cluding anti-particles, the complete information on the Lorentz structure of 
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the self-energy is available but in this case on has to cope with other problems 

uniiniiniiisi 

2 The Relativistic Brueckner Approach 
2.1 The Coupled Set of Equations 

In the relativistic Brueckner approach the nucleon inside the nuclear medium 
is viewed as a dressed particle in consequence of its two-body interaction 
with the surrounding nucleons. The in-medium interaction of the nucleons is 
treated in the ladder approximation of the relativistic Bethe-Salpeter (BS) 
equation 



where T denotes the T-matrix. V is the bare nuclcon-nucleon interaction. 
The intermediate off-shell nucleons in the scattering equation are described 
by a two-particle propagator iGG. The Pauli operator Q accounts for the 
influence of the medium by the Pauli-principle and projects the intermediate 
scattering states out of the Fermi sea. The Green's function G fulfills the 
Dyson equation 



Go denotes the free nucleon propagator while the influence of the surrounding 
nucleons is expressed by the nucleon self-energy E. In Brueckner theory this 
self-energy is determined by summing up the interaction with all the nucleons 
inside the Fermi sea in Hartree-Fock approximation 



The coupled set of equations ©-Q represents a self-consistency problem 
and has to be iterated until convergence is reached. 

Due to translational and rotational invariance, parity conservation and 
time reversal invariance the self-energy in isospin saturated nuclear matter 
has the general form S = Ss — "ff^E^. It depends on the Lorentz invariants 
k^, k- j and j^, with and being the baryon current and the nucleon four- 
momentum, respectively |18| . The invariants can also be expressed in terms 
of fco, |k;| and fcp, where k-p denotes the Fermi momentum. Furthermore the 
vector part of the self energy has contributions proportional to fc^ and to the 
current j'^. Defining the streaming velocity as = / \fp ^ the momentum 
k^ can be decomposed into contributions parallel and perpendicular to the 
streaming velocity, i.e. k^ ^ (k ■ u)u^ + A^'^k,y with the projector — 
g^" — u^u'^ . The vector part of the self-energy can then be written covariantly 

as iniini 




(2) 



G — Go + GqSG 



(3) 




So = —tr [j^u^S] = —tr [70 (9) 
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= SoU>' + K^'^'^K ■ (5) 
Thus the fuU self-energy reads 

r(fc, kp) = Ss{k, fcp) - [K{k, kp) uf' + U^ik, kp) A'^'^k,] (6) 
= Ss{k, fcp) - 70 Soik^kp) + 7 • kU^ik, fcp) |r_f (7) 

where the subscript RF indicates the respective expressions in the nuclear 
matter rest frame (u^ — 5'^°) |H1E|- The Ss, So and Sv components are 
Lorentz scalar functions which actually depend on fco,|k| and fcp- They follow 
from the self-energy matrix by taking the respective traces [llj 

Us = ^tr[S] (8) 

= ^^^Jr [zl-7.^. S] . ^tr [7 . k ^],, . (10) 

The Dirac equation for the in-medium spinor basis can be deduced from 
the Green's function. Written in terms of effective masses and momenta 

m* = M + ReS, , k*^ = k,, + Re Z"^ (11) 

the Dirac equation has the form 

[f-m* ~iImS]u{k)^0. (12) 

In the following we will work in the quasi-particle approximation and neglect 
the imaginary part of the self-energy from now on. Thus the effective nucleon 
four-momentum will be on mass shell even above the Fermi surface, fulfilling 
the relation k*^k*^ — m*^ . Since we only deal with the real part of the self- 
energy in the quasi-particle approximation we omit this in the notation. In 
the nuclear matter rest frame the four-momentum follows from Eq. (|ll|l 



k* = k(l + E^) , k*^E* = Vk2(l + rv)2+m*2 (13) 
which allows one to eliminate the Z'v-term in the Dirac equation, 

[{a-k)--f"m*]u{k)^E*u{k) , (14) 
by a rescaling of the effective mass and the kinetic energy 

E 



E* = — = v/k2 + m*2 . (15) 



1 + K 1 + 

The solution of the Dirac equation provide the in-medium nucleon spinor 
basis 
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Ux{k,kF) = d — , 2A|k| X\ , (16) 

V ^"^_F V B'(k)+r7i* / 



where £'*(k) = y/ls? + Th*p. xx denotes a two-component Pauli spinor with 
A = The normaUzation of the Dirac spinor is u\{k, kp)u\{k, kp) = 1. 
Since the in-medium spinor contains the reduced effective mass the matrix 
elements of the bare nucleon-nucleon interaction become density dependent. 
From the Dirac equation (|14() one derives the relativistic Hamiltonian, i.e. 
the single-particle potential U — j'^S. The expectation value of U, i.e. sand- 
wiching U between the effective spinor basis (fTT)|. yields the single particle 
potential 



which can be evaluated as 



kr.St" 



C/(fc,M = ^^.-^ (18) 



Vk2(l + rv)2+m*2 ^k^l + S^y 

In many applications 01201 single particle potential is only given in terms 
of a scalar and zero- vector component. This can be achieved by introducing 
reduced fields and So as 

y _ y A/f 

U, = m*-M= ; , K^E*-E = Uo-E*{k)U, .(20) 

The single particle potential has then the form 

U{k,kF) = ^Ss-So ■ (21) 

E* 

Frequently the reduced fields, Eq. I|2()(l . are used rather than the projected 
components since they represent the self-energy in a mean field or Hartree 
form. Thus they can easily be related to effective hadron mean field theory 
| 21l I22| . Such a representation is meaningful since the Z'v-contribution is a 
moderate correction. 



2.2 The In-Medium T-Matrix 

Before going into details I will shortly summarize the main assumptions which 
are made in the standard relativistic Brucckncr approach to solve the BS- 
equation 
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- No sea approximation. The subspace of negative energy states is omitted. 
In this way one avoids the dehcate problem of infinities which would gen- 
erally appear due to contributions from negative energy nucleons in the 
Dirac sea. The approximation is consistent with the usage of standard 
OBE potentials which are derived under the same assumption. 

- Thompson choice. The full two-body propagator iGG in the BS-equation 
is replaced by an effective two-body propagator propagator. The Thomp- 
son propagator (and similar the Blankenbecler sugar propagator) projects 
the intermediate nucleons onto positive energy states and restricts the ex- 
changed energy transfer by 5{k^) to zero. Thus the BS-equation is reduced 
to a three dimensional integral equation of the Lippmann-Schwinger type, 
the so called Thompson equation |2Hj . 

- Reference spectrums approximation. The momentum dependent effective 
mass rh* which enters into the Thompson propagator is replaced by an 
average value Trip (averaged over the Fermi sea). The approximation is 
justified as long as the self-energy exhibits a weak momentum dependence. 

- Angle-averaged Pauli operator. The Pauli operator is replaced by its angle- 
averaged counterpart which allows to solve the Thompson equation in a 
decoupled angular-momentum partial wave basis. 

- Quasi-particle approximation. The T-matrix is determined for on-shell 
scattering at the quasi-particle pole. Finite width spectral functions are 
not taken into account. 

In contrast to the self-energy, Eq. which has to be calculated in the 
nuclear matter rest frame, the Thompson Eq. H23|l is most naturally solved 
in the two-nucleon cm. frame. The Thompson propagator and similar the 
Blankenbecler-Sugar propagator imply that the time-like component of the 
momentum transfer in V and T is set equal to zero which is a natural con- 
straint in the cm. frame, however, not a covariant one. The Thompson equa- 
tion reads in the cm. frame 



where q = (qi — q2)/2 is the relative three-momentum of the initial state 
while k and p are the relative three-momenta of the intermediate and final 
states, respectively. The starting energy in Eq. is already fixed by -s/I* = 
2^*(q) = 2y/(i^ +rhf. If qi and q2 are nuclear matter rest frame momenta 
of the nucleons in the initial state, the boost-velocity u into the cm. frame 
is given by 



T(p,q,a;)| 



cm. 



^(P,q) 



(22) 
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u = P/VS* +P2 , (23) 

with the total three-momentum and the invariant mass P = qi + q2 and 
s* = (£'*(qi) + £'*(q2))^ - P^, respectively. In Eq. ^ x denotes the set of 
additional parameters x = {fcF,™J^, |u|} on which the T-matrix depends. 

The Pauli operator Q explicitely depends on the chosen frame, i.e., on the 
boost 3-velocity u into the cm. -frame. The Thompson equation (|23|l for the 
on-shell T-matrix (|p| — |q|) can be solved applying standard techniques de- 
scribed in detail by Erkelenz [21]. Doing so, one constructs the positive-energy 
helicity T-matrix elements from the IJAfLS*) -scheme. On-shell only five of the 
sixteen helicity matrix elements are independent which follows from general 
symmetries [S^. After a partial wave projection onto the | JMLS')-states the 
integral reduces to a one-dimensional integral over the relative momentum 
|k| and Eq. H23|l decouples into three subsystems of integral equations for 
the uncoupled spin singlet, the uncoupled spin triplet and the coupled triplet 
states. For this purpose the Pauli operator Q has to be replaced by an angle 
averaged Pauli operator Q fTSI . We are solving the integral equations by the 
matrix inversion techniques of Haftel and Tabakin |25j . Real and imaginary 
parts of the T-matrix are calculated separately by the principal- value treat- 
ment given in Ref. j26| . Due to the anti-symmetry of the two-fermion states 
the total isospin I of the two-nucleon system (1=0,1) can be restored by the 
selection rule: 

(_)L+s+i ^ (24) 

From the five independent on-shell amplitudes in the IJMLS*) -representation 
the five independent partial wave amplitudes in the helicity representation 
(for 1=0,1 and real and imaginary part separately) are obtained by inversion 
of Eq. (3.32) and then of Eq. (3.28) of Ref. [21. The summation over the 
total angular momentum J yields the full helicity matrix element 

diy{e){\Y>\\'^\'^\T'\x)\ |q|AiA2) = (pA'iA^Il3|T(a;)|qAiA2ll3). 

(25) 

Here 9 is the scattering angle between q and p and A = Ai — A2, A' = A'^ — Aj. 
The reduced rotation matrices d'{y{9) are those defined by Rose 27\. The 
matrix element H25(l is actually independent of the third component of the 
isospin I3. 

3 The Nuclear Self-Energy 

The easiest way to determine scalar and vector self-energy components di- 
rectly from the single particle potential U H17|l . Since U is obtained after 
complete summation over Dirac-indices of the T-matrix one by this way can 
avoid cumbersome projection techniques which are required using the trace 
formulas 0-1^01 ■ Using Eq. H21|l a fit to [/ delivers the density dependent 
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but momentum independent self-energy components and Sq. This method 
has e.g. been apphed in An attempt to extend this method and to ex- 
tract by fitting procedures momentum dependent fields '50] suffered by large 
uncertainties since one tries then to extract two functions out of one. 

A more accurate determination of the density and momentum dependence 
of the self-energy requires projection techniques for the in-medium T-matrix 
as outlined first by Horowitz and Serot (Ifl^ . That this procedure is also not 
free from ambiguities has been noticed relatively early |H1E|- The whole prob- 
lem arises from the no sea approximation in the standard approach. When 
calculations are performed in full Dirac space [lOL 1141 115L [T7] the Lorentz 
structure of the self-energy can uniquely determined from the information 
available form those matrix elements ((wjl^lu), (wlZ'lw)) which involve nega- 
tive energy states. The inclusion of negative energy excitations with 4 states 
for each spinor yields in total 4"^ = 256 two-body matrix elements for the 
T-matrix. Symmetry arguments reduce this to 44 for on-shell particles |28) . 
If one takes now only positive energy solutions into account this reduces to 
2* = 16 two-body matrix elements. 

For on-shell matrix elements the number of independent matrix elements 
can be further reduced by symmetry arguments down to 5. Thus, all on- 
shell two-body matrix elements can be expanded into five Lorentz invariants. 
But these five invariants are not uniquely determined since the Dirac matrices 
involve also negative energy states. The decomposition of a one-body operator 
into a Lorentz scalar and a Lorentz vector contributions depends therefore 
on the choice of these five Lorentz invariants. 

In nuclear matter the largest ambiguity arises concerning the determina- 
tion of pseudo-scalar {ps) and pseudo-vector [pv) T-matrix elements. The pv 
invariant in the medium is defined as 

with fcj, the initial and fc2 7 92 the final momenta of the scattering particles. 
For on-shell scattering of positive energy states the ps and pv matrix elements 
are identical (using the Dirac eq.) 

f4*~f\ 

u{q) Touip) = u{qh5u{p) . (27) 



. 2to* 

The ps vertex couples on the other hand maximally to negative energy 
whereas the pv vertex suppresses the coupling to antiparticles for on-shell 
scattering 

^(9) (^*) 75«(P) = . (28) 

To summarize: The ps and pv matrix elements of an on-shell two-body oper- 
ator, e.g. the T-matrix, can (in the positive energy sector) not be determined 
uniquely by projection techniques. If the matrix elements are known a priori^ 
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as e.g. for the bare NN interaction V, of course no problems arise. The same 
holds when the Dirac sea is included in the formalism. The full information 
on the T-matrix is available and the Lorentz structure of the self-energy is 
then uniquely determined |1(J[ I14L [TH ] . However, such an approach suffers 
from other problems (see next section). The ambiguity problem has e.g. been 








k [fm"'] 

Fig. 3. Influence of the various meson exchange contributions on the nucleon self- 
energy (scalar part). The solid line corresponds to the full calculation (Bonn A), 
the dotted line to auj and the dashed line to cruTr-exchange only. In all calculations 
the pv choice is used. 

pointed out by Nuppenau et al. ^Hl and ter Haar and Malfliet [S] proposed 
a recipe to cure which was used by various groups: perform the projection, 
take the ps matrix element and replace it (due to physical reasons) by a pv 
one. This procedure was called the pseudo-vector choice. Later on, this pro- 
cedure was critically examined by Fuchs et al. 1121 and it was shown that it 
completely fails in controlling the leading order 1 — 7r-exchange contribution. 

Fig. m shows the momentum dependence of the scalar self-energy compo- 
nent Ss at nuclear matter density p = 0.166 fm^^ obtained in the pv choice 
projection scheme, as it arises from the various meson exchange contribu- 
tions of the Bonn A potential. Taking only a and lo exchange into account 
the the momentum dependence is flat inside the Fermi sea. Including the 
pion we are already very close to the full DBHF result. The strong momen- 
tum dependence of the present calculation originates to a large extent from 
pion-exchange . 

It is a well known fact that a pseudo-scalar ttNN coupling leads to ex- 
tremely large pion contributions to the nuclear self-energy and contradicts 
soft pion theorems of ChPT . Therefore in OBEPs always V^r with pseudo- 
vector ttNN coupling is used. It is, however, instructive to test the pv choice 
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projection recipe for the case of the 7r-exchange JTJ- This is done in Fig. ^ 
where the Hartree-Fock self-energy from the l-7r-exchange potential (OPEP) 
is shown. Exact results for a ps and pv ttNN coupling can be compared to 
those obtained by projection techniques. It is seen that a description OPEP 
leads to extremely large self-energy components and a very strong momen- 
tum dependence. The pv coupling suppresses the contribution by nearly 
two orders of magnitude and even on that scale the momentum dependence 
is much less pronounced. At the Hartree-Fock level the K- self-energy can be 
computed directly or, alternatively, applying the same techniques as for the 
full T-matrix, i.e. going through the transformations from the \LSJ) basis to 
the helicity basis and finally via projection to the basis of covariant ampli- 
tudes. Doing so, it turns out that the pv choice projection fails to describe a 
pv pion exchange. Thus it is clear that the strong momentum dependence seen 
in the full self-energy (Fig. |2Jl is to large part due to spurious contributions 
form pseudo-scalar 7r-exchange. How these contributions can be eliminated 
has been discussed in |12| and ^31 and is explained in more detail below. 




Fig. 4. Hartree-Fock self-energy originating from OPEP. results for V-n with ps 
and pv coupling are compared to results obtained within the pv choice. Solid lines 
represent the scalar, dashed lines the vector self-energy. 



3.1 Covariant Representation of the T-Matrix 

To use the trace formulas, Eqs. H8I1U|I . one has to represent the T-matrix 
covariantly. A set of five linearly independent covariants is sufficient because 
on-shell only five helicity matrix elements appear as solution of the Thompson 
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equation. A linearly independent although not unique set of five covariants 
is given by the Dirac covariants 

S = 1 1, V 7^ ® 7^, T = ct'^'' ® a^^, A = 757^ ® 757^, P = 75 75- (29) 

Using this special set, dubbed as ps representation in the following, the on- 
shell T-matrix for definite isospin I can be represented covariantly as jj^l 

ri(|p|, 0, x) = Fl{\p\,e, x)s + Flj{\pie, x)y + F\M,e, x)t 

+ Fl(|p|,0,x)A + F]^(|p|,0,x)P . (30) 

Here p and 9 denote the relative three-momentum and the scattering angle 
between the scattered nucleons in the cm. frame, respectively. The direct 
(Hartree) amplitudes are are given by = and the exchange amplitudes 
(Fock) by = TT. The five covariant amplitudes can be obtained by matrix 
inversion of eq. (|30l) from the five helicity amplitudes T^. The nucleon self- 
energy in isospin saturated nuclear matter has then the form | 11| 

So.p{kM) = J (21^3 [AfIc.pFs+ r^pFy] , (31) 

where the isospin averaged amplitudes are defined as 

F,{\plO,x):=^[Ft''i\plO,x)+3Ft\\plO,x)] . (32) 

Eq. (|31|l shows that the self-energy can be expressed solely in terms of di- 
rect scalar and vector amplitudes Fsy, if these are derived from already 
anti-symmetrized helicity amplitudes (|25l) which obey the selection rule (I24II . 
Corresponding covariant exchange amplitudes F^{\p\,7T,x) are obtained from 
the exchange helicity amplitudes inverting a similar matrix for the exchange 
invariants S,V, T,A, P. The latter are related to the original invariants H29|l 
through a Fierz transformation. The same Fierz transformation relates also 
direct and exchange amplitudes -F]. Thus the direct scalar and vector ampli- 
tudes Fs^v in Eq. H31|l contain already contributions from all other exchange 
amplitudes pgllH] . 

Hence an explicit splitting of the already anti-symmetrized helicity ampli- 
tudes into direct and exchange parts i(T^(|p|, 0, x) — T^(|p|, tt, a;)) provides 
no additional information and is also not necessary. It becomes only relevant 
if one wants to replace the ps invariant by the pv invariant. However, now 
it becomes evident why the the above mentioned pv choice does not lead to 
the desired result: In this procedure the pseudo-scalar exchange amplitude 
i^p(|p|, TT, x) is kept fixed and interpreted as a pseudo-vector one, replacing 
the corresponding invariant P 1 — > PV. Since the other amplitudes, due to 
Fierz, contain also ps contributions such an replacement is incomplete. As 
can be seen from Fig.^the spurious ps contributions of the other amplitudes 
are still large. 
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To eliminate the such spurious pseudo-scalar contributions of the I-tt 
exchange completely one has to switch to another covariant representation 
of the T-matrix proposed by Tjon and Wallace [30] which we call full pv 
representation in the following 

T\\p\,9,x)^gl{\p\,e,x)S-gl{\p\,e,x)S+gi{\p\,9,x){A-A) 

+ gl.y{\p\,6,x)PV-g'~{\p\,e,x)PV . (33) 

The amplitudes g^{6) are explicitly given in p2|- In this scheme the pseudo- 
vector OPEP is exactly recovered. In Fig. [S] the corresponding self-energies 




Fig. 5. Range of uncertainty spanned by the various decomposition schemes of 
the T-matrix for: self-energy components (left); single particle potential (right). 
The nuclear matter density is chosen as p = 0.166 fm"'' and the Bonn A potential 
is used. 



obtained for the various decompositions are compared. Adopting the full pv 
representation, the space-like contribution turns out to be much smaller 
than in the ps or the standard pv choice. Therefore we show the reduced 
self-energies Eg and So in which is included for a better comparison. The 
pure ps and the full pv representation can be regarded as the limiting cases 
which give the range of uncertainty in the determination of the self-energy. 
The latter has the big advantage that this method ensures by construction 
a correct treatment of the pv OPEP at the Hartree-Fock level. Although 
the range of uncertainty is a few 100 MeV at the level of the self-energy 
components, it drops out to most extent for physical observables which are 
based on complete matrix elements where vector and scalar parts contribute 
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always with different sign. E.g., at the level of the single particle potential U 
the remaining uncertainty is only of about 10 20 MeV. 

3.2 Covariant Representations of the Subtracted T-Matrix 

The full pv representation successfully reproduces the HF nuclcon self-energy 
for the pion exchange with pv coupling. However, as pointed out in .12. . the 
full pv representation fails to reproduce the HF nucleon self-energy of other 
meson exchange potentials. Hence, it appears reasonable to treat the bare 
interaction V and the higher order ladder graphs of the T-matrix separately. 
Since the OBEPs are known analytically we can use a mixed representation 
of the form 

y = yj;,^ + K^l,,,, . (34) 

Here the tt- and 77-amplitudes are treated by the decomposition (|33|) while 
for the CT, w, /9, (5- amplitudes the representation (|30|l is applied. The higher 
order correlations of the T-matrix 

/oo „ 
VQGGT VWGGVr , (35) 

in the following called the subtracted T-matrix, can not be represented by 
such a mixed form since one can not disentangle the different meson con- 
tributions in the correlated ladder diagrams. The representation of the sub- 
tracted T-matrix remains therefore ambiguous. However, if the pion exchange 
contributes dominantly at the Hartree-Fock level a ps representation of the 
subtracted T-matrix should be more appropriate because of the higher order 
contributions of other meson exchanges. Thus the most favorable represen- 
tation of the T-matrix is given by the ps representation 

r^=Ti"«6+cr+<c.,p> • (36) 

Here the ps representation for T^^^ is determined via the matrix elements 

(pA;A;|Ti„,(a;)|qAiA2) {pXX\T\x) - V\x)\qXiX2) , (37) 

with subsequently applying the projection scheme as in Eq. An alter- 

native representation of the T-matrix is given by a representation 

t^''-tz+vi:;+vj:^^^, , (38) 

where the subtracted T-matrix is represented by the full pv representation 
1)33(1 . This representation is similar to the full pv representation of the full T- 
matrix, however, with the advantage that now the pseudo-scalar contributions 
in the bare nucleon-nucleon interaction, e.g. the l-co exchange potential, are 
represented correctly. 
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In ^21 the two representations and (|38|l for the higher order ladder 
graphs were studied in detail. These two representations set the range of the 
remaining ambiguity concerning the representation of the T-matrix, i.e. after 
separating the leading order contributions. The outcome is the following: 

- The dependence of the ladder kernel on the two different representation 
schemes is generally weak. This gives confidence that the ambiguities are 
to most extent removed as long as the leading Born term, in particular 
the pv OPEP, is treated correctly within the projection scheme. 

- The momentum dependence of the self-energy is moderate and close to the 
full pv case shown in Fig. [S] This observation also justifies the reference 
spectrums approximation. 

Larger differences between the ps and pv representations of the ladder kernel 
occur only at high densities. Here the ps representation appears to be more 
appropriate. Therefore we consider this as the best solution of the problem. 
More details can be found in In this context it should be noted that in 
| 31| the present subtraction scheme was extended to a special treatment of 
Vp where the tensor part leads also to a similar, but small on-shell ambiguity 
due to the Gordon identity. 

4 Nuclear Matter 

In this section I will discuss the role of correlations as well as implications 
of the various approximation schemes. Main emphasis will thereby be put on 
nuclear bulk properties, in particular on the nuclear equation-of-state (EOS), 
i.e. the binding energy per particle, and the nuclear saturation mechanism. 

4.1 The Equation-of-State 

In the relativistic Brueckner theory the energy per particle is defined as the 
kinetic plus half the potential energy 




In Fig. |Slthc EOSs, obtained in the various treatments, are compared. All 
calculations are based on the Bonn A interaction. First of all one sees that, 
except for a full ps treatment which is not correct for realistic potentials, 
the different calculations coincide at high densities. If one applies the pv 
choice the result is very close to that obtained by Brockmann and Machleidt 
(BM) [S]. This is somewhat surprising since there no projection scheme to 
the T-matrix has been applied but constant, i.e. momentum independent, 
self-energy components have been determined by a fit to the single parti- 
cle potential. As discussed in the previous sections the pv choice fails to 
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Fig. 6. Binding energy per particle as a function of nuclear matter density. As 
bare nucleon-nucleon interaction the Bonn A potential is used. For the T-matrix 
the subtraction scheme with the ps representation for the ladder kernel (solid) is 
compared to a ps representation of the full T-matrix (dash-dotted) , to the pv choice 
and to the result of |9j (BM, dotted). 

reproduce the pv OPEP contribution to the self-energy. One can estimate 
this effect at the level of the binding energy by the comparison of the ps 
representation for the subtracted and the ps representation for the full T- 
matrix. In the latter case the nucleons are less bound at small densities. The 
situation changes, however, around saturation density. The full ps represen- 
tation of the T-matrix contains maximal contributions from a pseudo-scalar 
ttNN coupling which leads to saturation properties closer to non-relativistic 
Brueckner calculations (see below). A correct pseudo- vector representation 
of the pion, as used in the subtraction scheme, suppresses this effect. Thus 
at smaller densities one obtains a larger binding, while around saturation 
density the EOS is more repulsive. 

In Fig. 13 we summarize the saturation points for iso-spin symmetric 
nuclear matter using different OBEPs as well as different approximation 
schemes. The saturation points for the two possible representations H36I38|I 
for the subtracted T-matrix are very similar . Therefore in the following I 
will consider the ps representation of the subtracted T-matrix as the optimal 
choice and compare this treatment with other works. With Bonn A one can 
reproduce the empirical saturation point of nuclear matter, shown as shaded 
region in the figure. The other Bonn potentials give less binding although 
the saturation density is always close to the empirically known value. Com- 
pared to the calculations of Brockmann and Machleidt [H] our Coester-line 
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is slightly shifted towards the empirical region which indicates that a refined 
treatment of the T-matrix leads to an enhancement of the binding energy 
connected with a reduced saturation density. In addition the result of ter 
Haar and Malfliet ^8 based on the Groningen OBEP is shown. All results 
were obtained in the no sea approximation. 
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Fig. 7. Saturation points of nuclear matter, obtained with different OBEPs and 
within different approximation schemes. Full symbols correspond to standard rela- 
tivistic Brueckner calculations, open symbols represent calculations which include 
Dirac sea contributions. The shaded area represents the empirical region of satura- 
tion. 



A second group of points (open symbols in Fig. [T]) includes explicitely 
negative energy states in the formalism. Here one has to keep in mind that 
standard OBEPs are derived within the no sea approximation and should 
therefore be used with caution in such calculations. The open triangle shows 
the result for Bonn C from which is close to the corresponding Bonn C 
values of the standard treatment. Significant differences occur when OBEPs 
are used which were derived in full Dirac space. This is the case for the cal- 
culations of Amorin and Tjon jj^ based on the Utrecht potential |22 and 
those of de Jong and Lenske ^H] based on the Gross OBEP "3^ (open tri- 
angles down). These results group at low densities which indicates strong 
additional repulsive components in the NN interaction resulting in a rather 
stiff EOS. As discussed in jT3] the coupling to the Dirac-sea, in particular 
the coupling to nucleon-antinucleon pairs (^-graphs) generates a strong dy- 
namical repulsion. In contrast to standard OBEPs, where repulsive Dirac 
sea contributions are effectively absorbed in a large lo coupling constant, 
OBEPs in full Dirac space generate such contributions dynamically. This is 
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also reflected in a significantly reduced to coupling constant. In nuclear mat- 
ter the Dirac sea contributions experience a medium dependence, primarily 
through the reduction of the mass gap. Consequently, the Z-graph contribu- 
tions are strongly enhanced at high densities which is the main source for the 
large repulsion observed in |14[ ITS] . Unfortunately, many-body calculations 
in full Dirac space show a strong sensitivity on off-shell effects, i.e. the corre- 
sponding form factors and to the three-dimensional reduction scheme of the 
BS-equation j^. Here certainly more efforts would be needed to control the 
influence of the Dirac sea in the many-body dynamics with higher accuracy. 
It should be noticed that in non-relativistic treatments Dirac sea contribu- 
tions can be accounted for on the level of three-body forces (see discussion 
below) . 

4.2 The Role of Correlations 

In order to examine the role of correlations it is instructive to compare the 
full DBHF theory to the mean field picture. In relativistic mean field theory 
(MFT) |18[I34| saturation occurs generally through the interplay between the 
large attractive scalar field Us, generated by the cr- meson, and the repulsive 
vector field So originating from the w-meson. In MFT the vector field grows 
linear with density while the scalar field saturates at large densities which 
leads finally to saturation. 
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Fig. 8. The effective nucleon mass in tlie DBHF approach, using different NN 
interactions and approximation schemes, is compared to mean field calculations. 
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As an illustrative example Fig. |S1 shows the density dependence of the 
effective nucleon mass m* for various models: We compare the present DBHF 
result (Bonn A, ps for the subtracted T-matrix) to the mean field result for 
Bonn A. In the latter only a and uj mesons contribute, exchange terms and 
contributions from other mesons vanish in iso-spin saturated nuclear matter 
^11^. In addition DBHF results from other groups are shown, i.e. those 
of in] (BM) , a more recent calculation of the Groningen group |SS| , and the 
original calculation of Horowitz and Serot (HS) ^HI using only a and uj 
exchange. The corresponding MFT result is shown as well. In the latter two 
calculations the coupling constants of QHD-I ^18. , adjusted to nuclear matter 
instead to NN scattering, have been used. Remarkable is that all calculations, 
though partially based on quite different models, lead to the same qualitative 
behavior. This indicates that the decrease of rh* and its tendency to saturate 
at high densities is dictated by relativistic dynamics. The different DBHF 
calculations lie thereby within a band of about 100 MeV which is set by the 
usage of different OBEPs and different approximation schemes. The reason 
for this common behavior of rh* is easy to understand: 

From eq. Ij^HI it follows that the scalar self-energy (jH)) is determined by 



In MFT the scalar amplitude Fg has to be replaced by the corresponding 
coupling constant of the scalar meson {gc/rricY. Eq. (|40|l represents then 
nothing else than the self-consistency equation for the effective mass 



which automatically leads to a saturating behavior for the attractive scalar 
field at large densities JS| . The momentum dependence of the T-matrix ele- 
ments is generally moderate This explains also why DBHF results can 
well be approximated within density dependent mean field theory |^ 
which means to replace Fs in by an average value Fs(A;f). Also the 
variation of such average amplitudes with density is in general moderate, 
however, with a tendency to decrease with density. Hence, one leading effect 
for saturation which takes place on the scale of the large scalar and vector 
fields of a few hundred MeV is present in full DBHF theory as well as in 
MFT. 

This does, however, not mean that the saturation mechanism is dominated 
by the mean field or Hartree contribution and exchange terms and higher 
correlations play only a minor role. The role of correlations can most easily 
be understood at the level of the two-nucleon wavefunction. As discussed e.g. 
in detail in |31 1^ and indicated in Fig. 13 correlations suppress the relative 
wavefunction at short distances. This reduces the short- and medium-range 
attraction of the (T-meson but even more effectively the short-range repulsion 
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Fig. 9. The EOS obtained for various interactions in the DBHF approach is com- 
pared to mean field theory (MFT). 



of the w-meson. Correlation lead therefore, first of all, to a general reduction 
of the magnitude of the self-energies and, secondly, to a significant reduction 
of the repulsive components of the interaction. As can be seen from Fig.|^ at 
saturation density the scalar field is reduced by about 100 MeV compared to 
MFT and even more at higher densities. The same holds for the vector field. 
For realistic OBEPs the quenching of the repulsive uj exchange is essential 
for the saturation mechanism. In a pure mean field picture the system turns 
e.g. out to be unbound for Bonn A. When couplings are already adjusted to 
nuclear matter in MFT, as done in QHD-I where (?5/47r — 10.84 is about 
half of the Bonn A value g'^/'in — 20 while is approximately the same, 
the higher order correlations lead to a significant softening of the EOS and 
shift the saturation point to higher densities [19^ . This behavior is illustrated 
in Fig. ini where DBHF results are compared to mean field calculations. The 
latter ones contain only contributions from a and uj exchange, however, with 
the coupling strengths of the corresponding OBEPs Bonn A and Groningen. 
The calculations denoted by (jlu in Fig. ^ are in both cases based on auj 
exchange only, however, now with the corresponding couplings of QHD-I. 



4.3 Role of the Pauli Operator 

Another important in-medium effect is represented by the Pauli operator Q 
which projects the intermediate states in the BS-equation ^ onto unoccu- 
pied phase space areas. The influence of Pauli blocking on the dynamics is 
most clearly seen if one considers directly matrix elements or, respectively in- 
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medium cross sections |371 EHl ED] ■ The differential on-shell cross section 
{p ~ q) is given by 

d<J=^^^\T{q,q,0)\'dn . (42) 

The squared matrix elements are obtained by the summation over 6 helicity 
helicity matrix elements (5 of them are independent) in the partial wave basis 
I^ElSni- From (|42|l one sees first of all that, compared to free scattering 
in the medium appears a suppression factor (to*/A/)^ which is solely due 
to kinematics. Furthermore, the Pauli operator modifies the optical theorem 
I4()j and damps in particular the imaginary part of the T- matrix which is 
directly proportional to Q (using R-matrix theory). For details see [10). 
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Fig. 10. Differential in-medium neutron-proton cross section for various densities 
at fixed laboratory energy of 250 MeV. The full results (thick lines) include partial 
waves up to J = 12 while the thin lines were obtained truncating higher partial 
waves. 

Another important Pauli effect is the suppression of higher partial waves 
in the two-body correlations. This effect can be estimated from Fig. IIUI where 
the medium dependence of the differential neutron-proton cross section (at 
■Blab = 250 MeV, using Bonn A) is displayed. The Fermi momenta correspond 
roughly to densities 1,2, and 3 (in terms of po = 0.1625 fm~^). The free np 
cross section is strongly forward-backward peaked. At moderate densities the 
presence of the medium tends to make the np differential cross section more 
isotropic. At backward angles the cross section are decreasing with density. 
At forward angles the behavior is more complicated: At moderate densities 
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the cross section is reduced but at high densities {p = 2 and 3 po) a strong 
enhancement of the forward scattering amphtude can be observed. Similar 
results have been obtained by Li and Machleidt At 3po the cross section 
turns out to be again highly anisotropic and to be dominated by a p-wave 
component. There occurs generally a suppression of higher partial waves with 
increasing density: At po one needs partial waves up to at least J < 6 to 
approximate the full result (J = 12), at 2po the partial waves J < 4 are 
almost sufficient and at Spo the behavior is dominated by s + p-waves. 

It is quite natural that Pauli blocking is most efficient for the low momen- 
tum components of the interaction, generated mainly by iterated 7r-exchange, 
while high momentum components from iterated heavy meson exchange 
(tr, Ld) arc much less affected. 

Another effect which is closely related to the Pauli operator is a possible 
onset of superfluidity at low densities. E.g. in the finite temperature approach 
of Aim et al. [35] a critical enhancement of the np cross section at low densities 
has been observed which was attributed to the onset of superfluidity. Crucial 
for such a superfluid state are contributions from hole-hole scattering in the 
Pauli operator [Q — (1 — / — /)) which are absent in the standard Brueckner 
approach (Q = (1 — /)(1 — /)). However, as discussed in 0J a signature of a 
bound pair state can appear at low densities even when hole-hole scattering 
is neglected in the Pauli operator. In ^Hl such an resonance like enhancement 
was seen in the amplitudes which correspond to the quantum numbers of the 
deuteron, i.e. the ^Si, ^Di and the ^Si-^Di transition channels. Therefore 
the low density enhancement of the np cross section can be interpreted as a 
precursor of a superfluid state. The same effect has been discussed in ^T] . 

5 Relativistic versus Non-Relativistic BHF 

In contrast to relativistic DBHF calculations which came up in the late 80ies 
non-relativistic BHF theory has already almost half a century's history. The 
first numerical calculations for nuclear matter were carried out by Brueckner 
and Gammel in 1958 '5. Despite strong efforts invested in the development 
of methods to solve the Bethe-Goldstonc (BG) equation, the non-relativistic 
counterpart of the BS equation, it turned out that, although such calcula- 
tions were able to describe the nuclear saturation mechanism qualitatively, 
they failed quantitatively. Systematic studies for a large variety of NN inter- 
actions showed that saturation points were always allocated on a so-called 
Coester-line in the E/A — p plane which does not meet the empirical region 
of saturation. In particular modern OBEPs lead to strong over-binding and 
too large saturation densities where relativistic calculations do a much better 
job. Several reasons have been discussed in the literature in order to explain 
the success of the relativistic treatment. In the following I will recapitulate 
the main arguments for this difference. 
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5.1 Continuous Choice versus Gap Choice 

Brueckner theory converges in terms of the hole-line expansion (for a recent 
review see 0). In lowest order Brueckner theory (2 hole- lines) the effective 
2-particle propagator in the BG-equation leads to a gap at the Fermi surface 



Gl2 — 



e(q) - Eik) 



(43) 



with e(q) = U{kF) + the single particle energies below the Fermi mo- 
mentum (starting energy) and E(k) — the energy of the intermediate 
states above fcp- From (|23() it is evident that the relativistic propagator does 
not contain such a gap. The continuous choice advocated by the Liege group 
|H] assumes the single particle potential to be valid also above fcp. This is 
in line with the relativistic propagator where fields are present below and 
above the Fermi momentum. Compared to the gap (or standard) choice the 
continuous choice shifts the Coester-line significantly towards the empirical 
region 021 E3|- It was further shown by the Catania group jjlj that at the 
3-hole-line level both choices lead to almost identical results. In the continu- 
ous choice already lowest order Brueckner theory (2-hole-lines) is very close 
to the result of the 3-hole-line expansion which suggests a faster convergence 
of the continuous choice. 
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Fig. 11. Saturation points of relativistic (full symbols) versus non-relativistic (open 
symbols) BHF calculations (continuous choice). 
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5.2 Dirac Effects and Quenching of the Tensor Forces 

The saturation mechanisms in relativistic and non-relativistic theories are 
quite different. In relativistic MFT the vector field grows linear with den- 
sity while the scalar field saturates at large densities. In MFT this is the 
essential effect which leads to saturation. As discussed above, the density 
dependence of the scalar and vector DBHF self-energy is similar to MFT, 
however, exchange contributions and correlations lead a significant reduction 
of their absolute magnitude. A pure mean field picture works when a and lo 
couplings are adjusted to nuclear matter, however, when fitted to free scat- 
tering data no saturation may occur at the mean field level, depending on 
the choice of the interaction. Exchange contributions and correlations are of 
crucial importance in order to obtain the saturation point at a scale which 
is set by the binding energy. A genuine feature of relativity, present in MFT 
as well as in DBHF, is the presence of large scalar and vector fields of a few 
hundred MeV size and a strongly decreasing effective nucleon mass with the 
tendency to saturate at high densities. 

In the non-relativistic case the situation is quite different: The saturation 
mechanism takes place exclusively on the scale of the binding energy, i.e. a 
few ten MeV. It cannot be understood by the absence of a tensor force [32] . In 
particular the second order OPEP is large and attractive at high densities and 
its interplay with Pauli-blocking leads finally to saturation. Relativistically 
the interaction Y receives a density dependence since matrix elements of 
V are built between in-medium spinors H16|l . As discussed in [T6) the tensor 
interactions 14 , and the Pauli part of Vp experience an additional reduction 
by the factor {fh* /MY. This reduction is at a first glance not completely 
obvious since it enters through the momentum dependence of second and 
higher order pseudo- vector OPE (and similar for 77 and p) where the range of 
the intermediate momentum k is controlled by ifi* . Consequently, the OPEP 
contribution to the binding energy is significantly reduced in the relativistic 
approach. Consistent with this observation is the present result (Figs. El and 

obtained by a pure pseudo-scalar representation of the T-matrix within 
the projection scheme. As discussed above this leads to stronger weights of the 
self-energy contributions from pion exchange and in turn to stronger binding 
at higher density. It should, however, be kept in mind that in this calculation 
only the representation of the T-matrix is taken as ps, the 14- itself is still 
used with pv coupling. 

In summary we have two genuine features of relativistic dynamics which 
are closely connected, act, however, at different scales: the saturation of the 
scalar attraction takes place on the scale of the large self-energy fields and 
the quenching of the tensor force at the scale of the binding energy. As argued 
in 4fia 42, the latter effect is probably responsible for the improved Coester 
lines compared to BHF. In the language of effective field theory it may be 
tempting to relate this two scales with chiral fluctuations on top of large 
background fields originating from QCD condensates |47j . Remarkably, such 
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an approach can finally lead to very similar self-energies as the present DBHF 
calculations |47 |. 

5.3 Resonance Degrees of Preedom and Three-Body Forces 

Since the inclusion of explicit resonance degrees of freedom (DoFs) in the 
formalism is closely related to the occurrence of 3-body forces (3-BFs) I will 
discuss here both aspects in combination. 

The most important resonance is of course the Z\(1232) isobar. At low 
and intermediate energies it provides large part of the intermediate range 
attraction and generates most of the inelasticity above the pion threshold. 
Intermediate A states appear in elastic NN scattering only in combination 
with at least two-isovector-meson exchange (tttt, Trp, . . .) and give rise to a 
new class of box diagrams. As has been shown by the Bonn group this class 
of diagrams can satisfactorily be absorbed into the effective cr-exchange. If the 
A is maintained as an explicit DoF in TV TV scattering, it provides additional 
attraction and the corresponding a strength has to be readjusted. This leads 
e.g. in the work of the Groningen group to a reduction of g^/^i^ = 7.4 (w/o 
A) to gl/Ai: = 6.4 including A DoFs. 

In many-body calculations explicit A DoFs give rise to additional satura- 
tion, shifting the saturation point away from the empirical region. This can 
be understood in the following way while the elementary aNN vertex 
is not modified in the medium, corresponding box diagrams with intermedi- 
ate A and nucleon lines are affected. Dressing of the propagators and Pauli 
blocking of the nucleon state suppresses their contribution. Iterated to all or- 
ders, the maintenance of explicit A DoFs (instead of a stronger a exchange) 
results therefore in less attraction. Quantitatively this effect has in detail 
been studied by ter Haar and Malfliet |H]. Fig. ll2lshows the result of a more 
refined calculation from the Groningen group ^Hj which includes also the A 
self-energy. As one sees, the loss of binding energy is quite substantial. As 
pointed out in Refs. the inclusion of non-nucleonic degrees of freedom 

has to be performed with caution: Freezing out resonance DoFs generates 
automatically a class of three-body forces which contains nucleon-resonance 
excitations. 3-BFs with intermediate A excitations provide again a strong in- 
termediate range attraction, A^*(1440) excitations lead to small net repulsion 
|3T]. Hence there exist strong cancellation effects between the repulsion due 
to box diagrams and contributions from 3-BFs. Such ring type diagrams of 
third and forth order in the hole-line expansion have e.g. been conducted in 
|52| . A consistent treatment requires therefore to consider non-nucleonic DoFs 
and many-body forces on the same footing ^H] • Many-body forces which are 
exclusively based on nucleon degrees of freedom can systematically be gener- 
ated within chiral perturbation theory. Next to leading order all 3-BFs cancel 
[55] while non- vanishing contributions appear at NNLO |54| . 

Another class of 3-BFs which has extensively been studied within non- 
relativistic BHF involves virtual excitations of nucleon-antinucleon pairs. 
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Fig. 12. Saturation points of DBHF calculations from the Groningen group without 
and including explicit A degrees of freedom 49 , and of BHF calculations (50| based 
on the AVis potential without and including 3-body forces. 



Such Z-graphs are in net repulsive and lead to a shift of the satura- 
tion point away from the non-relativistic Coester line towards its relativistic 
counterpart where the DBHF results are allocated. The calculation shown in 
Fig. 1121 includes both, NN and as well as nucleon-resonance excitations pT) . 

It is often argued that in non-relativistic treatments 3-BFs play in some 
sense an equivalent role as the dressing of the two-body interaction by in- 
medium spinors in Dirac phenomenology. Both mechanisms lead indeed to 
an effective density dependent two-body interaction V which is, however, of 
different origin. In the medium 3-BFs can be considered as a renormaliza- 
tion of the meson vertices and propagators. Z-graphs are explicitely included 
when DBHF calculations are performed in full Dirac space, in the no sea 
approximation they are in some way effectively included through the usage 
of OBEPs with large ui couplings. 



6 Summary 

An overview on the present status of relativistic Brueckner calculations for the 
nuclear many-body problem was given. Using modern one-boson-exchange 
potentials such calculations provide a qualitatively satisfying - and param- 
eter free - description of the nuclear saturation mechanism. Concerning the 
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extractions of the precise nuclear self-energy and its Lorentz structure there 
arise on-shell ambiguities due to lack of information on two-body matrix 
elements in full Dirac space when the approach is restricted to the positive 
energy sector (no sea approximation). A method to minimize the correspond- 
ing uncertainties was discussed. 

Similar to relativistic mean field theory, Dirac phenomenology together 
with the structure of the NN interaction extracted from free scattering im- 
plies the existence of large scalar and vector fields. However, exchange con- 
tributions and higher order correlations reduce the magnitude of these fields 
compared to MFT and are essential for a quantitative saturation mechanism. 
When calculations are performed in full Dirac space part of the repulsion is 
generated from sea excitations which requires to renormalize the NN poten- 
tials. However, in the latter case the many-body dynamics can presently not 
be controlled with the same accuracy as in the standard approach based on 
the no sea approximation. In non-relativistic BHF nucleon-antinucleon ex- 
citations can be accounted on the level of three-body forces which leads to 
qualitatively similar results as in relativistic approach with solely two-body 
interactions. However, a consistent treatment of 3-BFs is a subtle problem 
which is closely connected to the introduction of non-nucleonic degrees of 
freedom, i.e. nuclear resonances. A future perspective would be the applica- 
tion of chiral NN potentials [SS] where 3-BFs can consistently be by power 
counting. Another somewhat complementary challenge, relevant for the ap- 
plication to relativistic heavy ion reactions, is an extension to higher energies, 
using thereby high precision potentials above the pion threshold |5f)| . 
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